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We  report  on  the  measurements  of  the  shock  equation  of  state  (Hugoniot)  of  an  Al/Fe203/epoxy 
composite,  prepared  by  epoxy  cast  curing  of  powder  mixtures.  Explosive  loading,  with  Baratol, 
trinitrotoluene  (TNT),  and  Octol,  was  used  for  performing  experiments  at  higher  pressures,  in  which 
case  shock  velocities  were  measured  in  the  samples  and  aluminum,  copper,  or  polymethyl 
methacrylate  (PMMA)  donor  material,  using  piezoelectric  pins.  The  explosive  loading  of  the  metal 
donors  (aluminum  and  copper)  will  be  discussed.  Gas  gun  experiments  provide  complementary 
lower  pressure  data  in  which  piezoelectric  polyvinylidene  fluoride  (PVDF)  stress  gauges  were  used 
to  measure  the  input  and  propagated  stress  wave  prohles  in  the  sample  and  the  corresponding  shock 
propagation  velocity.  The  results  of  the  Hugoniot  equation  of  state  are  compared  with  mesoscale 
hnite-element  simulations,  which  show  good  agreement.  ©  2007  American  Institute  of  Physics. 
[DOI:  10.1063/1.2719272] 

I.  INTRODUCTION 

Equation  of  state  studies  on  epoxy-based  particulate 
composite  systems,  particularly  alumina-epoxy  systems, 
have  been  reported  in  the  literature,  as  epoxy-based  com¬ 
posites  are  becoming  more  widely  used  for  structural  appli¬ 
cations  because  of  their  low  densities  and  relatively  high 
strengths.  The  properties  of  the  composite  can  be  tuned 
based  on  the  properties  and  volume  fraction  of  the  particulate 
phase(s).  Results  for  shock  wave  propagation  in  these  com¬ 
posite  materials  show  dispersive  waves  with  no  elastic  pre- 
cursor.  In  particular,  Anderson,  et  al.  and  Deas  et  al.  have 
observed  that  the  Hugoniot  data  for  these  composite  systems 
showed  a  linear  relationship  between  pressure  and  particle 
velocity  over  the  range  of  conditions  investigated. 

The  research  presented  in  this  article  focuses  on  equation 
of  state  measurements  performed  on  Al  and  Fe203  powders 
cast-cured  in  an  epoxy  binder.  The  combination  of  aluminum 
and  iron  oxide  has  been  studied  as  the  traditional  thermite 
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reaction.  The  shock  compression  behavior  of  these  powder 
mixtures  was  hrst  reported  by  Holman  et  a/.^  While  no  evi¬ 
dence  of  chemical  reaction  was  revealed  in  these  powder 
mixtures  at  pressures  less  than  15  GPa,  it  was  found  that  the 
shock  compression  response  was  dominated  by  the  dissimilar 
compressibilities  of  the  softer  Al  and  harder  Fe203  powders. 

In  this  article,  equation  of  state  measurements  are  pre¬ 
sented  for  aluminum  and  iron  oxide  powders  in  an  epoxy 
matrix.  Results  over  a  wide  range  of  pressures  are  obtained 
from  both  gas  gun  and  explosive  loading  experiments.  A 
discussion  of  the  explosive  loading  of  metal  plates  used  as 
donor  materials  is  also  included.  Einally,  the  experimental 
Hugoniot  is  compared  with  that  calculated  by  mesoscale 
hnite-element  simulation.^ 

II.  EXPERIMENTAL  PROCEDURE 
A.  Sample  preparation 

Spherical  aluminum  powder  and  blocky  iron  oxide  pow¬ 
der  were  obtained  in  particle  sizes  of  approximately  2  /u,m 
for  the  Al  and  0.3  /u-m  for  Pe203,  as  shown  in  Figs.  1(a)  and 
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FIG.  1.  Micrographs  of  (a)  spherical  aluminum  powder  (1-2  /rm  particle  size);  (b)  blocky  iron  oxide  powder  ~0.3  /rm,  and  (c)  as-cast  composite  material 
showing  Al,  Fe203  agglomerates  (~  8  /rm  average  size)  in  an  epoxy  matrix. 


1(b),  respectively.  The  powders  were  blended  with  Epon  826 
epoxy  (Hexion  Specialty  Chemicals)  cured  with  diethanola¬ 
mine  (DEA)  hardener,  poured  into  aluminum  molds,  vacuum 
degassed  at  70  °C,  and  cast  in  the  form  of  50  mm  diameter 
by  125  mm  long  rods.  Sections  of  ~5  mm  thickness  were 
sliced  from  the  rod  and  used  for  performing  Hugoniot  mea¬ 
surements.  The  weight  percents  of  the  Al,  Fe203,  and  epoxy 
were  12.63,  37.37,  and  50%,  respectively.  Prior  to  perform¬ 
ing  the  Hugoniot  measurements  the  samples  were  character¬ 
ized  to  determine  their  density,  acoustic  wave  speed,  and 
microstructure.  Figure  1(c)  shows  the  microstructure  of  a 
polished  surface  of  the  sample,  revealing  the  rounded  Al  par¬ 
ticles  and  Fe203  agglomerates  (of  ~8  /zm  average  size)  in¬ 
terspersed  in  an  epoxy  matrix. 

B.  Gas  gun  loading  experiments 

Four  equation  of  state  experiments,  listed  in  Table  I, 
were  conducted  at  Georgia  Tech  using  the  80  mm  diameter 
single-stage  light-gas  gun  at  impact  velocities  between  500 
and  950  m/s.  The  experiments  were  instrumented  with  Bauer 
piezoelectric  polyvinylidene  fluoride  (PVDF)  stress  gauges® 
obtained  from  Ktech  Corporation,  Albuquerque,  NM.  The 
gauges  were  mounted  on  the  front  and  rear  surfaces  of  the 
sample  (as  shown  in  Fig.  2),  which  was  sandwiched  between 
a  copper  driver  and  a  PMMA  backer  plate.  An  aluminum 
projectile  with  a  copper  flyer  plate  (9.5  mm  nominal  thick¬ 
ness),  backed  by  an  air  gap,  was  used  to  impact  a  copper 
driver  plate  (6.3  mm  nominal  thickness).  The  gauge  mounted 
between  the  copper  driver  and  the  sample  provided  an  “in¬ 
put”  stress  profile,  and  the  gauge  mounted  between  the 
sample  and  the  PMMA  backer  provided  the  “propagated” 
stress  profile.  The  thicknesses  of  the  flyer,  driver,  and  target 
plates  were  adjusted  to  ensure  that  a  steady-state  shock  wave 
propagated  through  the  sample  thickness  without  attenuation. 
The  travel  time  between  the  two  gauges  (less  the  time  of 
travel  through  the  gauge  package)  was  used  to  determine  the 


shock  wave  speed  in  a  sample  of  known  thickness.  Assuming 
hydrodynamic  conditions  and  steady-state  wave  propagation 
through  the  sample,  the  measured  input  stress  and  the  shock 
wave  speed,  combined  with  the  jump  conditions,  were  used 
to  determine  the  particle  velocity  and  pressure  Hugoniot. 

C.  Explosive  loading  experiments 

Twenty-five  experiments,  detailed  in  Table  II  and  shown 
schematically  in  Fig.  3,  were  conducted  at  Eglin  AFB  using 
explosive  plane  wave  lenses  (PWL)  to  generate  shock  load¬ 
ing  at  higher  pressures.  In  these  experiments,  the  PWL  was 
in  contact  with  an  explosive  pad  of  Baratol,  TNT,  or  Octol 
(12.5  or  25.4  mm  thick)  to  generate  a  range  of  pressures  in  a 
6061-T6  aluminum  (Al),  OFHC  copper  (Cu),  or  PMMA  do¬ 
nor  plate.  In  experiments  JJH8-12,  the  explosive  pad  was 
12.5  mm  thick,  and  in  the  other  experiments  the  explosive 
pad  was  25.4  mm  thick. 

In  explosive  loading  experiments,  the  shock  wave  travels 
from  the  donor  to  the  sample.  The  shock  velocity  in  the  Al  or 
Cu  donor  material  and  the  sample  was  measured  using  pi¬ 
ezoelectric  pins  (Dynasen  CA-1135),  which  were  placed  in 
holes  drilled  into  the  sample  at  differing  depths.  Using  the 
shock  velocity  in  the  donor  plate  and  the  sample,  the  remain¬ 
ing  Hugoniot  properties  for  the  sample  were  determined  us¬ 
ing  impedance  matching. 

In  experiments  JJH8-12,  six  piezoelectric  pins  were 
placed  in  the  sample  and  the  donor  plate  to  characterize  the 
shock  velocity  in  both  materials.  In  experiments  JJH 13-26, 
nine  piezoelectric  pins  were  used.  Experiments  JJH36-52 
were  designed  to  characterize  the  explosive-metal  interaction 
[Fig.  3(b)].  In  these  experiments,  20  piezoelectric  pins  were 
used  to  get  a  statistical  representation  of  the  shock  velocity 
in  the  metal  donor  in  contact  with  the  explosive. 

A  representative  plane  wave  lens  similar  to  those  used  in 
these  experiments,  was  characterized  to  determine  the  pla¬ 
narity  of  the  lens  in  the  confined  configuration  of  the  experi- 


TABLE  I.  Experimental  details  for  equation  of  state  experiments  utilizing  gas  gun  loading. 


Expt. 

no. 

Target  thickness 
(mm) 

Target  density 
(g/cm^) 

Cu  flyer  thickness 
(mm) 

Impact  velocity 
(m/s) 

Copper  donor 
thickness  (mm) 

0303 

3.9131  ±0.0053 

1.7900  ±0.0 122 

9.5 123  ±0.0025 

523±27 

6.2814±0.0025 

0308 

3.8034±0.0066 

1.8645±0.0019 

9.4996  ±0.0025 

553±43 

6.3017±0.0025 

0311 

3.8212±0.0038 

1.8503±0.0082 

9.5204  ±0.0025 

714±89 

6.3251  ±0.0051 

0403 

3.8735±0.0064 

1.8588±0.1765 

9.5148±  0.0025 

944  ±4 

6.3119±0.0051 
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Gas  gun  barrel  (80  mm  dia)  Yfilocity  pins 


PMMA  backer 
(9  mm) 


Trigger 

pins 

PMMA  window 
(17.5  mm) 

Backer  PVDF 
gauge 


Epoxy-cured 
Al-bFe203  sample 
(3.8  mm) 


FIG.  2.  Schematic  of  experiments  performed  on  Al/Fe203/epoxy  compos¬ 
ites  using  gas  gun  loading  showing  the  placement  of  stress  gauges  used  to 
monitor  input  and  propagated  stress  profiles. 


merits.  The  streak  camera  record  showing  the  break¬ 
out  of  the  detonation  wave  is  illustrated  in  the  inset  of  Fig.  4. 
The  breakout  of  the  wave  was  digitized  as  a  function  of 
position  and  time,  as  shown  in  Fig.  4;  the  simultaneity  of  the 
plane  wave  lens  breakout  is  ±40  ns. 


III.  RESULTS  AND  DISCUSSION 

In  all  explosive  loading  experiments,  the  shock  velocity 
was  calculated  by  fitting  a  line  to  the  position  versus  time 
data  for  all  pins.  The  shock  wave  velocity  {U^  for  a  given 
sample  are  provided  in  Tables  III  and  IV.  A  tilt  correction 
was  applied  to  the  shock  velocity  data  based  on  a  method 
described  by  Dick.^  In  this  method,  the  pin  arrival  times,  r,-. 


TABLE  II.  Experimental  details  for  equation  of  state  experiments  utilizing 
explosive  plane  wave  lens  loading. 


Expt. 

no. 

Explosive 

Thickness 

(mm) 

Donor 

Sample 

JJH8 

TNT 

12.5 

Al 

Al/Fe203/epoxy 

JJH9 

TNT 

12.5 

PMMA/Al 

Al/Fe203/epoxy 

JJHIO 

TNT 

12.5 

PMMA 

Al/Fe203/epoxy 

JJHll 

Octol 

12.5 

Al 

Al/Fe203/epoxy 

JJH12 

Octol 

12.5 

PMMA/  Al 

Al/Fe203/epoxy 

JJH13 

Octol 

12.5 

PMMA 

Al/Fe203  /  epoxy 

JJH16 

Baratol 

25.4 

Al 

Donor  only 

JJH42 

Baratol 

25.4 

Al 

Donor  only 

JJH43 

Baratol 

25.4 

Al 

Donor  only 

JJH25 

Baratol 

25.4 

Cu 

Al  /  Fe2 0 3  /  epoxy 

JJH51 

Baratol 

25.4 

Cu 

Donor  only 

JJH52 

Baratol 

25.4 

Cu 

Donor  only 

JJH14 

TNT 

25.4 

Al 

Donor  only 

JJH22 

TNT 

25.4 

Al 

Donor  only 

JJH36 

TNT 

25.4 

Al 

Donor  only 

JJH37 

TNT 

25.4 

Al 

Donor  only 

JJH24 

TNT 

25.4 

Cu 

Al/Fe203/epoxy 

JJH45 

TNT 

25.4 

Cu 

Donor  only 

JJH46 

TNT 

25.4 

Cu 

Donor  only 

JJH15 

Octol 

25.4 

Al 

Donor  only 

JJH39 

Octol 

25.4 

Al 

Donor  only 

JJH40 

Octol 

25.4 

Al 

Donor  only 

JJH26 

Octol 

25.4 

Cu 

Al/Fe203/epoxy 

JJH48 

Octol 

25.4 

Cu 

Donor  only 

JJH49 

Octol 

25.4 

Cu 

Donor  only 

Piezoelectric  Pin$ 


Piezoelectric  Pins 


LEGEND 

06061-  T6  A1 

^  TNT,  Octol  or  Baratol  Pad 

0TNT 

0PMMA 

USample  of  Interest 

[^Steel 

[^Donor.6061.T6  Al, 

On!CCii,orPMMA 

r~n  Octol 


FIG.  3.  (Color  online)  Schematic  of  equation  of  state  experiments  on 
Al/Fe203/epoxy  composites  using  explosive  loading  where  (a)  is  the  con¬ 
figuration  used  for  experiments  JJH8-16  and  24-26  and  (b)  is  the  configu¬ 
ration  used  for  JJH36-52. 


the  pin  depths,  z,,  and  the  pin  coordinates  on  the  surface  of 
the  donor  plate,  x,  andy^,  are  fitted  to  the  equation 

tj  =  Pi  + P2Zi  + PiXi  + Fiji,  (1) 

where  the  subscript  i  indicates  each  pin  time  and  correspond¬ 
ing  pin  location.  The  constants  P3  and  P4  are  a  measure  of 
the  tilt  of  the  shock  wave  with  respect  to  the  plane  of  the  pin 
circle.  The  time  the  pins  would  have  discharged  had  there 
been  no  tilt,  t^j,  can  be  represented  as 

hi  =  h  -  PiXi  -  P^yi  =  Pi  +  P2Zi,  (2) 

where  Pj  is  the  intercept  and  P2  is  the  slope  of  the  corrected 
time,  t^.  The  reciprocal  of  this  slope  is  the  corrected  shock 
velocity,  recorded  in  Tables  III  and  IV  as  “11^  Tilt  Corr.”  In 
experiments  with  only  one  pin  circle  in  the  donor  plate  and 
one  in  the  sample,  the  correction  was  performed  to  the  center 
of  the  pin  circle.  In  experiments  (JJH39-52)  with  two  pin 
circles,  the  correction  was  performed  to  the  center  of  the 
explosive  charge. 


® 

E 


u> 

3 

E 


Position  (mm) 


FIG.  4.  Breakout  of  detonation  wave  from  TNT/Octol  plane  wave  lens  as  a 
function  of  position  and  time.  The  inset  picture  shows  the  original  streak 
camera  record. 
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TABLE  III.  Measured  shock  velocity  (U^  with  calculated  pressure  (P)  and  particle  velocity  (1/^)  in  aluminum 
and  copper  donors  in  contact  with  TNT,  Octol,  and  Baratol  compared  with  literature  values.  In  JJH8-12,  the 
explosive  pad  was  12.5  mm  thick.  In  all  other  experiments,  the  explosive  pad  was  25.4  mm  thick. 


Expt. 

Explosive 

Donor 

V, 

u. 

Us 

Us 

P 

Up 

no. 

(mm//Lts) 

tilt  corr. 

avg. 

lit. 

(GPa) 

(mm//,ts) 

(mm/^ts) 

(mm//xs) 

(mm/^ts) 

IJHIO 

TNT 

PMMA 

5.42±0.19 

5.42±0.19 

12.0±0.9 

L86±0.13 

JJH9 

TNT 

PMMA  /  A1 

6.97±0.52 

6.99±0.48 

23.2±6.9 

L23±0.36 

JJH8 

TNT 

A1 

6.62  ±0.44 

6.63±0.42 

17.1±5.8 

0.95  ±0.32 

JJH13 

Octol 

PMMA 

5.72±0.25 

5.75±0.15 

14.2±0.8 

2.08±0.10 

JJH12 

Octol 

PMMA/  A1 

7.29±0.29 

7.33±0.15 

29.2±2.3 

L48±0.11 

IIHll 

Octol 

A1 

7.12±0.12 

7.13±0.08 

25.5  ±1.2 

L32±0.06 

JJH16 

Baratol 

A1 

6.73±0.67 

6.99±0.40 

6.7 

JJH42 

Baratol 

A1 

6.29  ±0.40 

6.74±0.06 

2L2±2.2 

L14±0.10 

JJH43 

Baratol 

A1 

6.33±0.45 

6.90±0.02 

6.89±0.20 

JJH25 

Baratol 

Cu 

4.48±0.52 

4.79±0.20 

4.79±0.09 

JJH51 

Baratol 

Cu 

4.32±0.31 

4.71  ±0.07 

4.5 

24.5±2.9 

0.57±0.06 

JJH52 

Baratol 

Cu 

4.55±0.29 

4.88±0.05 

JJH14 

TNT 

A1 

6.30±0.49 

6.39±0.39 

JJH22 

TNT 

A1 

6.93±0.80 

7.31  ±0.46 

6.89±0.20 

7.0 

2L5±3.0 

L15±0.16 

JJH36 

TNT 

A1 

6.60±0.35 

6.90±0.11 

JJH37 

TNT 

A1 

6.67±0.33 

6.95±0.06 

JJH24 

TNT 

Cu 

4.93±0.12 

4.95±0.03 

JJH45 

TNT 

Cu 

4.56±0.19 

4.69±0.05 

4.79±0.14 

4.8 

24.5±4.8 

0.57±0.09 

JJH46 

TNT 

Cu 

4.57±0.21 

4.73±0.01 

JJH15 

Octol 

A1 

6.66±0.64 

6.81±0.50 

7.8 

JJH39 

Octol 

A1 

7.38±0.20 

7.45±0.12 

L38±0.25 

JJH40 

Octol 

A1 

7.24±0.19 

7.32±0.06 

7.19±0.34 

26.8±6.0 

JJH26 

Octol 

Cu 

4.64±0.40 

4.66±0.38 

JJH48 

Octol 

Cu 

5.13±0.11 

5.16±0.02 

4.91  ±0.25 

5.4 

28.9±8.9 

0.65±0.17 

JJH49 

Octol 

Cu 

4.87±0.10 

4.91±0.01 

In  determining  the  shock  velocity  (for  both  uncorrected 
and  corrected  values)  the  displacement  data  vector  (z^ata)  is  a 
linear  function  of  time,  i.e.,  z(t)  =  ait-i-a2>  where  and  <32 


are  unknown  coefficients.  A  straightforward  maximum  like¬ 
lihood  estimator  (MLE)  approach  is  implemented^  where  the 
vector  of  error  in  the  data  (i.e.,  the  residual). 


TABLE  IV.  Shock  velocity  (I/,),  pressure  (P),  and  particle  velocity  {Uj)  in  Al/Ee203/epoxy  samples  from  gas 
gun  and  explosive  loading.  The  experimentally  measured  values  are  italic. 


Expt. 

no. 

Initial  density 
(g/cm^) 

Us 

(mm/^ts) 

Us 

tilt  COIT. 

(mm/ /xs) 

Up 

(mm//xs) 

P-m 

(GPa) 

^OUt 

(GPa) 

0303 

L790±0.012 

3. 17  ±0.01 

0.372±0.006 

2.11  ±0.02 

1.76±0.22 

0308 

L865±0.002 

3.47±0.01 

0.364±0.004 

2.35±0.35 

0311 

L850±0.008 

3.65±0.01 

0.499  ±0.004 

3.37±0.I2 

2.95±0.06 

0403 

L856±0.177 

3.92±0.01 

0.641  ±0.005 

4.66±0.31 

3.56±0.4 

JJH25 

L850±0.027 

3.42±0.26 

3.50±0.15 

0.97±0.24 

6.3±L6 

JJH26 

L850±0.027 

4.84±0.13 

4.85  ±0.06 

0.78±0.41 

7.0±3.6 

JJH24 

L850±0.027 

4.18±0.09 

4. 19  ±0.04 

L13±0.02 

8.7±0.3 

JJH8 

L850±0.027 

5.16±0.11 

5.16±0.11 

L22±0.41 

1L7±3.9 

JJH9 

1.850±0mi 

4.7 1  ±0.36 

4.77±0.26 

L63±0.49 

14.4±4.4 

IJHIO 

L850±0.027 

5.12±0.23 

5.14±0.18 

L57±0.11 

15.0±L2 

JJHll 

L850±0.027 

5.34±0.16 

5.35±0.13 

L71±0.08 

16.9±0.9 

JJH13 

L850±0.027 

5.46±0.46 

5.46±0.46 

L76±0.12 

17.7±L9 

JJH12 

L850±0.027 

4.93±0.33 

4.96±0.27 

L96±0.16 

17.9±L8 

JJH22/23 

L850±0.027 

5.47±0.36 

5.52±0.24 

1L88±0.46 

19.2±4.8 
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r  =  z(a,r)  -Zdata, 


(3) 


is  minimized  in  a  least-squares  sense  to  find  the  best  estimate 
of  the  unknowns  (x),  i.e., 

a  ^  argmin^(r  ^n;V),  (4) 

where  11^  is  the  covariance  matrix  of  the  uncertainties  in 
the  z(f)  data.^  However,  is  not  known  prior  to  the  analy¬ 
sis,  so  the  uncertainty  in  the  data  of  length  is  estimated 
from  an  unweighted  solution  of  Eq.  (4).**’  Here  11^=1  ,  such 
that 


r^r 


(5) 


Up 


and  correspondingly 

n,  =  a^I.  (6) 

The  estimate  of  the  uncertainty  in  x  is  then  found  by  apply¬ 
ing  a  QR  decomposition  to  the  Jacobian  matrix,  J,  for  the 
displacement  function  with  respect  to  the  estimation 
parameters, 

J  =  [t  1]  =  QR,  (7) 

where  Q  is  unitary  and  R  is  upper  triangular.  The  estimated 
covariance  of  the  parameters  at  the  best  estimate  is  then 
found  from 

n,=  (R%R)-'.  (8) 

Using  this  approach,  the  error  and  correlations  in  the  shock 
velocity,  both  corrected  and  uncorrected  can  be  determined. 

For  the  donor  materials,  the  remaining  parameters  (par¬ 
ticle  velocity,  C/„,  and  pressure,  P)  can  be  determined  based 
on  the  published  Hugoniot  data,  ’  which  are  given  by 

U,  =  5.35+l.34Up,  (9) 

for  6061-T6  Al,  with  an  initial  density,  po=8.93  g/cm^,  and 

U,=  3.94+lA9Up,  (10) 

for  OFHC  copper,  with  an  initial  density,  po=8.93  g/cm^, 

and 


U,  =  2.6Q+l.52Up,  (11) 

for  PMMA,  with  an  initial  density,  po=1.19  g/cm^,  and  the 

12 

Rankine-Hugoniot  jump  conditions,  namely. 


FIG.  5.  Schematic  of  impedance  matching  method  where  line  1  is  the  Hugo¬ 
niot  for  the  donor  material,  line  2  is  the  reflected  Hugoniot  for  the  donor 
material,  line  3  has  the  slope  PoUl  for  the  sample  material,  and  point  4  is  the 
shock  state  in  the  sample  material. 


In  the  sample  material,  the  shock  velocity  is  determined 
using  the  same  method  as  described  for  the  donor  material. 
To  determine  the  remaining  parameters,  the  impedance 
matching  technique,  shown  graphically  in  Fig.  5,  was  used. 
The  particle  velocity  in  the  sample  is  determined  from  the 
intersection  of  the  reflected  donor  Hugoniot  and  a  line  with  a 
slope  of  PqUI  for  the  sample.  The  particle  velocity  in  the 
sample.  Up,  can  be  determined  from  the  following  equation, 
in  which  the  superscript  s  indicates  the  sample  and  d  indi¬ 
cates  the  donor: 


,  -  b  ±  'Jb^  -  4ac 

K  = - > 

‘  2a 

(15a) 

where 

a  =  PqS‘', 

(15b) 

b=-{4pis‘‘u‘;,+pict+p^uD, 

(15c) 

c  =  4p^5''(U^)2  +  2p^C^t/^. 

(15d) 

The  particle  velocity  could  also  be  approximated  by  estimat¬ 
ing  the  reflected  Hugoniot  for  the  donor  as  a  straight  line 
with  a  slope  of  poU^-  Using  this  estimation,  the  error  associ¬ 
ated  with  the  sample  particle  velocity  can  be  determined 
from 


P  =  PoU,Up. 


(12) 


Equations  (9)-(ll)  take  the  form  Us=CQ+SUp.  Here,  Cq  rep¬ 
resents  the  bulk  sound  speed  and  the  S  represents  a  linear 
empirical  dependence  of  the  shock  velocity  on  the  particle 
velocity.  The  errors  in  the  particle  velocity  and  pressure  are 
propagated  according  to  standard  quadrature  estimates. 


AUp 


AU.  \ 

U,-Co/’ 


(13) 


and 


AP  =  F 


H- 


(14) 


AUi=  \[^AUty  + 


[\ip^ou:+ptut)  r 


+  \-^AUlf, 


(16a) 


where 

/  2pipiu:ui  \ 

\(p-^t/j  +  p^uf)V 


(16b) 


The  pressure  in  the  sample  and  its  associated  error  can  be 
determined  from  Eqs.  (12)  and  (14),  in  which  the  density, 
shock  velocity,  and  particle  velocity  are  those  associated 
with  the  sample  material  rather  than  the  donor  material. 
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A.  Explosive-donor  interaction 

Using  the  combination  of  Baratol,  TNT,  and  Octol  pads 
with  copper  and  aluminum  driver  plates,  shock  velocities 
ranging  from  5-7  km/s  can  be  achieved  in  the  metal  donor 
material.  This  provides  a  range  of  input  pressures  into  the 
samples  and  enables  a  reasonable  quantification  of  the  Hugo- 
niot  in  the  mid  to  high  pressure  range. 

Table  III  presents  the  measured  shock  velocity  and  cal¬ 
culated  pressure  and  particle  velocity  for  all  explosive-donor 
combinations.  Experiments  JJH8-13  were  conducted  with  a 
12.7  mm  thick  explosive  pad  in  contact  with  the  plane  wave 
lens.  It  is  believed  that  this  is  not  a  large  enough  distance  to 
achieve  a  steady  detonation  wave  in  the  explosive  pad,  so 
later  experiments  utilized  an  explosive  pad  that  was  25.4  mm 
thick.  These  data  remain  usable,  although  the  pressure  is 
likely  below  that  of  the  predicted  detonation  pressure  in  the 
explosives. 

It  can  be  seen  that  the  experiments  are  very  repeatable. 
However,  the  earlier  experiments,  JJH14-16  and  JJH24-26 
do  not  agree  as  well  with  the  better-calibrated  experiments 
that  were  performed  later.  This  could  be  due  to  the  increased 
number  of  pins  in  later  experiments  (20  vs  9)  and  the  im¬ 
proved  tilt  correction  and  statistics  associated  with  the  in¬ 
creased  number  of  pins.  Additionally,  the  pins  in  the  early 
experiments  used  a  thicker  piezoelectric  crystal  (0.5  vs  0.25 
mm),  which  have  a  longer  fall-off  time.  The  pins  are  placed 
in  close  proximity  so  the  shock  velocity  is  measured  over  a 
small  distance.  In  the  early  experiments,  several  pin  records 
overlapped,  making  it  difficult  to  distinguish  the  time  when 
the  signal  deviates  from  the  baseline,  and,  subsequently,  add¬ 
ing  to  the  uncertainty  in  the  experiment. 

1.  Equation  of  state  measurements 

The  measured  results  (shaded)  and  calculated  values  ob¬ 
tained  from  the  gas  gun  and  explosive  loading  experiments 
are  presented  in  Table  IV.  The  shock  velocity  values  corre¬ 
spond  to  those  measured  based  on  time  of  travel  through  the 
sample  thickness  between  the  input  and  propagated  stress 
gauges  in  the  gas  gun  experiments  and  between  piezoelectric 
pins  in  explosive  loading  experiments.  Stress  profiles  re¬ 
corded  by  the  PVDF  gauges  in  the  gas  gun  experiments  pro¬ 
vided  the  magnitude  of  the  stress  entering  the  sample.  Figure 
6(a)  shows  examples  of  stress  profiles  recorded  by  the  input 
stress  gauge  as  the  shock  wave  generated  upon  impact  of  the 
flyer  on  to  the  driver  enters  the  sample.  Figure  6(b)  shows 
the  transmitted  stress  pulses  for  experiments  0303,  0311,  and 
0403  following  travel  through  approximately  4  mm  thickness 
of  composite  material,  with  corresponding  peak  pressures  of 
1.7,  2.9,  and  3.6  GPa,  respectively.  The  rise  time  to  peak 
pressure  decreases  with  increasing  stress.  The  lower  pressure 
experiment  (0303)  shows  a  highly  dispersed  stress  wave  pro¬ 
file  due  to  wave  dispersion,  similar  to  behavior  observed  by 
Anderson  et  al.  in  alumina-filled  epoxy.  The  material  tested 
in  experiment  0303  had  a  lower  density  than  the  other 
samples  tested.  This  point  is  included  for  comparison,  but  is 
not  used  when  analyzing  the  data. 

The  results  combining  the  data  obtained  from  the  gas 
gun  and  explosive  loading  experiments  are  plotted  in  the 


FIG.  6.  (Color  online)  PVDF  gauge  traces  from  (a)  input  gauge  and  (b) 
transmitted  gauge. 


U^-Up  and  P-Up  space,  as  shown  in  Figs.  7(a)  and  7(b).  The 
line  fitted  through  the  data  points  (except  for  the  datum  cor¬ 
responding  to  shot  0303,  which  had  a  sample  of  lower  den¬ 
sity)  exhibits  a  smooth  trend  without  any  obvious  disconti¬ 
nuities. 

The  error  in  the  explosive  loading  experiments,  detailed 
in  Table  IV,  is  larger  than  the  —3%  error  in  stress  and  0.2% 
error  in  timing  measurements  in  the  case  of  the  gas  gun 
experiments.  These  larger  errors  could  be  due  to  the  curva¬ 
ture  of  the  shock  wave  after  it  enters  the  donor.  Additionally, 
the  0.5  mm  crystal  piezoelectric  pins  were  used  in  all  of 
these  experiments,  and  deconvolution  of  the  individual  pin 
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(a)  Up  (km/s) 


(b)  Up  (km/s) 

FIG.  7.  (Color  online)  Hugoniot  plots  based  on  calculated  and  measured 
values  obtained  from  gas  gun  and  explosive  loading  experiments:  (a)  pres¬ 
sure  vs  particle  velocity;  (b)  shock  velocity  vs  particle  velocity,  in  compari¬ 
son  with  mesoscale  hnite-element  calculations.  The  linear  fit  to  this  plot 
dehnes  the  dynamic  bulk  sound  speed  (3.18  km/s)  and  the  empirical  param¬ 
eter  S  (1.16). 


responses  proved  to  be  difficult  in  several  cases.  Another 
source  of  error,  especially  for  data  obtained  using  a  copper 
donor,  is  the  large  impedance  mismatch  between  the  copper 
and  composite  material,  which  may  translate  to  larger  uncer¬ 
tainty  compared  to  the  6061-T6  Al.  In  both  the  gas  gun  and 
explosive  loading  experiments,  the  errors  are  associated  with 
measurements  of  the  experimental  parameters  and  are  not 
meant  to  represent  multiple  tests  at  the  same  input  pressure. 
In  the  gas  gun  experiments,  these  errors  are  due  to  the  inher¬ 
ent  limitations  of  the  measurement  devices,  and,  in  the  ex¬ 


plosive  loading  experiments,  these  errors  are  due  to  the  lin¬ 
ear  fit  to  the  time-distance  data  required  for  determining  the 
shock  velocity. 

Considering  a  linear  ht  to  the  data  shown  in  the  shock 
velocity  versus  particle  velocity  plot  in  Fig.  7(a),  values  of 
the  dynamic  bulk  sound  speed,  Cq,  and  empirical  parameter, 
S,  were  obtained  for  the  Al/Fe203/ epoxy  cast  cured 
samples.  For  these  materials,  values  of  Cq  equal  to  3.17  km/s 
and  S  equal  to  1.16  were  obtained.  The  bulk  sound  speed  of 
the  material  from  ultrasonic  measurements  is  2.2  km/s,'^  and 
the  values  calculated  from  mass  averaging  the  properties  of 
the  components  are  a  Cq  value  equal  to  4.41  km/s  and  S 
equal  to  1.38.*^  The  curve-fit  depicted  in  Fig.  7(b)  is  the 
translation  of  the  linear  relationship  determined  in  the  Us- Up 
space  to  P-Up  space.  This  assumes  that  the  stress  measured 
by  the  PVDF  gauges  is  equivalent  to  the  hydrodynamic  re¬ 
sponse,  which  was  shown  by  Deas  et  al?  to  be  true  for 
alumina-epoxy  composites  tested  in  a  flyer  impact  configu¬ 
ration  with  velocities  ranging  from  190-670  m/s. 

Although  there  is  considerable  scatter  in  the  data,  par¬ 
ticularly  in  the  intermediate  particle  velocity  range,  the 
smooth  trend  exhibited  by  the  data  from  both  sets  of  experi¬ 
ments,  with  no  apparent  discontinuities  in  the  Hugoniot,  in¬ 
dicates  that  there  are  no  chemical  or  physical  changes  occur¬ 
ring  at  pressures  up  to  19  GPa  in  the  Al/Fe203/epoxy 
composite  samples  containing  the  large  volume  fraction  of 
epoxy. 

2.  Comparison  of  measured  Hugoniot  with 
finite-etement  modei  simutations 

Shock  wave  propagation  in  discrete  particles  mixtures 
has  been  studied  at  the  mesoscale  by  Benson  and  co-workers 
for  a  range  of  material  systems,  e.g.,  Cu  powders, 
HMX,^®  Ti-SiC,*^  and  Nb-Si.^*’'^  Recently,  mesoscopic 
finite-element  models  have  been  developed  to  enable  the 
study  of  shock  waves  in  Al/Fe203/epoxy  composite 
systems.  In  the  aforementioned  modeling  efforts,  thermo¬ 
mechanical  responses  have  been  resolved  at  the  level  of  the 
microstructure.  The  following  discussion  provides  a  broad 
overview  of  the  methodology  for  modeling  shock  wave 
propagation  in  Al/Fe203/ epoxy  composite  systems;  the  lit¬ 
erature  should  be  consulted  for  in-depth  modeling  details. 

The  shock  wave  simulations  are  performed  in  a  2D  mul¬ 
timaterial  Eulerian  finite-element  code  developed  by 
Benson.  Eulerian  formulations  are  necessary  in  these  cal¬ 
culations  due  to  the  severity  of  deformation  (material  flow). 
The  governing  equations  (conservation  of  mass,  momentum, 
and  energy)  are  solved  sequentially  using  an  operator  split.^^ 
Here,  the  central  difference  method  is  used  to  advance  the 
solution  in  time  during  the  Lagrangian  step.  Advection  algo¬ 
rithms  are  used  to  transport  material  between  cells  of  the 
spatially  fixed  Eulerian  mesh  during  the  Eulerian  step  (solu¬ 
tion  remap)  while  time  remains  constant.  Heat  generation 
and  heat  transfer  are  included  so  that  thermal  and  mechanical 
effects  are  coupled  in  the  solution.  The  plane-strain  condition 
is  imposed  due  to  the  computational  demands  of  fully  3D 
calculations. 

The  explicit  microstructure  of  the  particle  mixtures  must 
be  reconstructed  for  simulation.  Discrete  sets  of  particles  are 
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generated  from  lognormal  size  distributions  that  correspond 
to  those  measured  in  actual  particle  mixtures.  Volume  frac¬ 
tions  of  the  particle  phases  and  pores  in  the  simulations  are 
matched  to  those  studied  in  experiments.  Particles  are  lo¬ 
cated  in  the  domain  using  a  random  sequential  addition 
process.  A  two-point  spatial  correlation  is  imposed  on  the 
distribution  of  A1  particles;  here,  a  simulated  annealing  tech¬ 
nique  is  used  to  evolve  the  first-nearest-neighbor  distribu¬ 
tions  in  the  A1  phase  to  that  estimated  in  the  actual  micro¬ 
structure.  The  final  result  is  a  2D  rectangular  cross  section  of 
a  particle  mixture  to  be  analyzed  using  finite-element  simu¬ 
lation. 

Boundary  conditions  are  applied  to  the  domain  to 
achieve  nominal  ID  shock  wave  propagation.  A  velocity 
boundary  condition  is  imposed  on  one  edge  of  the  domain  to 
generate  the  shock  wave.  The  boundary  conditions  for  all 
remaining  edges  are  “rollers”  (i.e.,  normal  displacements  are 
constrained  to  be  zero  and  tangential  displacements  are  un¬ 
constrained).  The  simulations  are  terminated  before  the 
shock  front  reaches  the  remote  boundary  so  that  wave  reflec¬ 
tions  are  avoided. 

A  set  of  constitutive  models  is  required  to  define  the 
stress-strain  behavior  of  each  phase.  During  simulation,  the 
stress  response  is  decomposed  into  hydrostatic  and  deviatoric 
terms,  which  necessitates  the  definition  of  an  equation  of 
state  (EOS)  and  strength  model  for  each  phase.  The  Mie- 
Gruneisen  EOS  is  used  to  model  the  hydrostatic  response  of 
the  A1  and  epoxy  phases,  while  the  Birch  Murnaghan  EOS 
(Ref.  24)  is  used  for  the  hydrostatic  response  of  the  Ee203 
phase.  The  Klepaczko  model^^  defines  the  rate-dependent 
plastic  flow  and  strain  hardening  of  the  A1  phase.  The  Hasan- 
Boyce  model^^  defines  the  rate-dependent  behavior  of  the 
epoxy  matrix.  Since  a  physically  based  constitutive  stress- 
strain  model  is  not  available  for  the  Ee203  phase,  a  simple 
isotropic  elastic -plastic  strength  model  with  a  small  degree 
of  work  hardening  is  adopted. 

Results  are  extracted  from  the  simulations  by  using 
simple  homogenization  techniques  to  calculate  EOS  data  for 
the  macroscopic  mixture.  The  shock  velocity  in  a  simulation 
is  calculated  by  tracking  the  position  of  the  shock  front  with 
respect  to  time.  The  shock  pressure  is  calculated  by  volume 
averaging  the  pressure  in  a  region  of  the  microstructure  that 
is  sufficiently  far  away  from  both  the  applied  velocity  bound¬ 
ary  condition  and  the  shock  wave  front.  The  particle  velocity 
is  known  in  the  simulations;  it  is  imposed  (as  a  boundary 
condition)  to  match  that  measured  in  an  experiment. 

The  U^-  Up  and  P-  Up  relations  calculated  from  the  model 
are  plotted  in  Eigs.  7(a)  and  7(b),  respectively.  As  shown  in 
Eig.  7(a),  shock  wave  velocity  calculations  agree  well  with 
experimental  data  in  the  high-velocity  regime  {Up 
>1.000  km/s),  as  differences  are  less  than  8%;  differences 
are  somewhat  larger  (capped  at  14%)  in  the  low- velocity 
regime  (t/p<  1.000  km/s).  This  trend  is  in  agreement  with 
prior  observations  by  Benson'^  for  Cu  powders,  where  larger 
deviations  in  the  shock  wave  speed  calculations  were  ob¬ 
served  at  lower  particle  velocities.  As  shown  in  Eig.  7(b), 
stationary  pressure  calculations  show  good  agreement  with 
experimental  data  in  the  low-velocity  regime,  as  differences 
are  less  than  7%;  differences  are  somewhat  larger  (capped  at 


12%)  in  the  high-velocity  regime.  Despite  discrepancies  in 
the  U^-Up  relation,  it  is  encouraging  to  observe  a  linear  de¬ 
pendence  between  the  shock  wave  velocity  and  particle  ve¬ 
locity,  as  such  relations  describe  many  material  systems.  The 
deviations  in  calculations  of  t/j  and  P^t  are  tentatively  attrib¬ 
uted  to  the  2D  approximation  of  the  micro  structure. 

IV.  CONCLUSIONS 

Equation  of  state  experiments  in  the  2-19  GPa  pressure 
range  were  conducted  on  Al/Ee203/epoxy  composites. 
Lower  pressure  experiments  were  performed  using  an  80  mm 
diameter  single-stage  light-gas  gun.  The  input  and  propa¬ 
gated  pressures,  as  well  as  the  shock  velocity  in  the  material 
were  recorded  using  PVDE  piezoelectric  gauges.  In  order  to 
achieve  higher  pressures,  explosive  loading  of  the  samples 
was  conducted  using  TNT/Octol  plane  wave  lenses  with 
Baratol,  TNT,  and  Octol  pads  and  aluminum,  copper,  and 
PMMA  donor  plates  to  vary  the  input  pressure.  Shock  veloc¬ 
ity  measurements  were  obtained  using  piezoelectric  pins  in 
the  sample  and  in  donor  materials,  in  order  to  perform  im¬ 
pedance  matching  to  determine  the  remaining  properties  in 
the  sample.  The  data  from  both  sets  of  experiments  illustrate 
a  continuous  trend  with  no  evidence  of  any  physical  or 
chemical  change  occurring  in  the  Al/Pe203/epoxy  compos¬ 
ites  at  pressures  up  to  19  GPa.  It  is  satisfying  that  the  gas  gun 
data  and  the  explosively  derived  data  are  comparable.  Ei- 
nally,  mesoscale  finite-element  modeling  results  show  good 
agreement  with  experimental  data. 
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